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ABSTRACT 

Shock waves driven by the release of energy at the center of a cold ideal gas sphere of initial density 
p cx approach a self-similar behavior, with velocity R cx R s , as R —> oo. For u> > 3 the solutions 
are of the second-type, i.e., S is determined by the requirement that the flow should include a sonic 
point. No solution satisfying this requirement exists, however, in the 3 < uj < uj g (^) "gap" (uj g = 3.26 
for adiabatic index 7 = 5/3). We argue that second- type solutions should not be required in general 
to include a sonic point. Rather, it is sufficient to require the existence of a characteristic line r c (t), 
such that the energy in the region r c (t) < r < R approaches a constant as R — > 00, and an asymptotic 
solution given by the self-similar solution at r c {t) < r < R and deviating from it at r < r c may 
be constructed. The two requirements coincide for u) > uj g an d the latter iden tifies S = solutions 
as the asymptotic solutions for 3 < to < uj g (as suggested by iGruzinovl 120031 ) . In these solutions, 
r c is a Co characteristic. It is difficult to check, using numerical solutions of the hydrodynamic 
equations, whether the flow indeed approaches a S — self-similar behavior as R — » 00, due to the 
slow convergence to self-similarity for uj ~ 3. We show that in this case the flow may be described by 
a modified self-similar solution, dla R/dhi R = S with slowly varying S(R), rj = d5/dhiR <C 1, and 
spatial profiles given by a sum of the self-similar solution corresponding to the instantaneous value 
of 5 and a self-similar correction linear in r\. The modified self-similar solutions provide an excellent 
approximation to numerical solutions obtained for uj ~ 3 at large R, with 5 —> (and rj ^ 0) for 

3 < UJ < LJg. 

Subject headings: hydrodynamics — self-similar — shock waves 



1. INTRODUCTION 

Self-similar solutions to the hydrodynamic equations 
describing adiabatic one-dimensional flows of an ideal gas 
are of interest for mainly two reasons. First, the nonlin- 
ear partial differential hydrodynamic equations are re- 
duced for self-similar flows to ordinary differential equa- 
tions, which greatly simplifies the mathematical problem 
of solving the equations and in certain cases allows one 
to find analytic solutions. Second, self-similar solutions 
often describe the limiting behavior approached asymp- 
totically by flows which take place over a c haracteristic 
scale, R, which diverges o r tends to ze ro Jsee lSedovlUQSl 
IZel'dovich fc Raizedfl967l: iBaren blatt 1991, for reviews) . 

It is reasonable to assume that in the limit R — > 
oo(0) the flow becomes independent of any character- 
istic length scale. Using dimensional arguments, it is 
possible to show that in thi s case the flow fields must 
be of the self-similar for m (jZel'dovich fc Raizcr 1967; 
IWaxman fc Shvartsll20i0l ) 

u(r, f) - R£U($), c(r, t) = R£C(0, p{r, t) = BR e G(£), 

where u, c, and p are the fluid velocity, sound speed, and 
density, respectively (the pressure is given by p = pc 2 /'*/), 
and 

R = AR S , £(r,t) = r/R(t) (2) 

(for a somewhat different appro ach to self-similarity, 
based on Lie group met h ods, see [Co ggcshal l fc Axford 
119861: ICoggeshalll 119911: ICoggeshall fc Mever-ter-Vehn 
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I1992D . For a self-similar solution of the form given 
in Equatio ns ([! ]) and ([2]), the hydrodynamic equations, 
Equations (|A1|) ■ are replaced with a single ordinary dif- 
ferential equation, Equation (IA2[) . 

dU _ Ai(t/,C) 
dC ~ A 2 (C/,<7)' 

and one quadrature, Equation (|A3I) . 

d\n£, _ A(U,C) dln£ _ A(U,C) 

dU ~ Ai (U,C) ° r dC ~ A 2 (U,C)' 

A, Ai, and A 2 are given by Equa tions (IA6I) . As 
illustrated in Section [2 1 (see also Guderley 1942; 
| Mever-ter-Vehn fc Schalkl |1982t IWaxman fc Shvartsl 
Il993f ). many of the properties of self-similar flows may 
be inferred by analyzing the contours in the (U, C)-planc 
determined by Equation (IA2|) . 

In this paper, we revisit the "strong explosion prob- 
lem" , which is one of the most familiar problems where 
asymptotic self-similarity is encountered. Consider the 
blast wave produced by the deposition of energy E within 
a region of characteristic size d at the center of an ini- 
tially cold (p = at r > d) gas sphere with initial density 
po = Kr~ M (at r > d). As the shock radius R diverges, 
we expect the flow to approach a self-similar solution of 
the form given by Equations ([]} and ([2]). The asymptotic 
flow is describ ed by the Sedov-von Neumann-Taylor 
(ST) solutions (|Sedovl 119461; Ivon Neumann] 119471; iTavlorl 
|1950() for uj < 3, an d by the solutions derived by W ax- 
man & Shvrats (WS: IWaxman fc Shvartslll993l l2010h for 
uj > 3. The ST solutions describe decelerating shocks 
(S = (uj — 3)/2 < 0) and are of the "first-type", where 




u c (4/3) ~ 8.22 
w c (5/3) ~ 7.69 



w,(5/3) ~ 3.26 
w„(4/3) ~ 3.13 
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Fig. 1. — Self-similar exponent 8 as a function of lu for 7 = 4/3 
(dashed line) and 7 = 5/3 (solid line). 



the similarity exponents, 8 and e, are determined by di- 
mensional considerations. The WS solutions describe ac- 
celerating shocks (8 > 0) and are of the "second-type" , 
where the similarity exponents are determined by the 
condition that the solutions must pass through a singu- 
lar point of Equation (|A2|) . 

We revisit the strong explosion problem for several 
reasons. First, there exists a "gap" in the (7, w)-plane, 
3 < co < ujgij), where neither the ST nor the WS solu- 
tions describe the asymptotic flow (u> g is increasing with 
7, uig = 3 for 7 = 1 and co a ~ 3.2 6 for adiabatic index 
7 = 5/3; Wax man fc Shvartslll993h . Our first goal is to 
close this "gap". We argue in Section [5] that second- type 
solutions should not be required in general to include a 
sonic point, and that it is sufficient to require the ex- 
istence of a characteristic line r c (t), such that the en- 
ergy in the region r c (t) < r < R approaches a constant 
as R — > 00. We show that the two requirements coin- 
cide for ui > co g and that the latter requirement iden- 
tifies 8 = solutions as the asymptotic solutions for 
3 < co < co a . T his result is in agreement with that of 
iGruzinovl (|2003l ). who suggested based on heuristic argu- 
ments that the R oc t solutions are the correct asymptotic 
solutions in the gap. As we explain in some detail at the 
end of Scction [2~2l the validity of the heuristic arguments 
is not obvious. We use a dif ferent reasoning, bas e d on a n 
extension of the analysis of I Waxman fe Shvartsl (|1993f) . 

In Section [3~T1 we compare the asymptotic, R/d 3> 1, 
behavior of numerical solutions of the hydrodynamic 
equations, Equations (|A1|) . to that expected based on 
the 8 = self-similar solutions. We show that the con- 
vergence to self-similarity is very slow for co ~ 3. Hence, 
it is difficult to check using numerical solutions whether 
the flow indeed approaches a 8 — self-similar behav- 
ior as R — > 00. We show in Section 13.21 that in this 
case the flow may be described by a modified self-similar 
solution, d In R/d In R = 8 with slowly varying 8(R), 
r] = d8/d\nR <C 1, and spatial profiles given by a sum 
of the self-similar solution corresponding to the instanta- 
neous value of 8 and a self-similar correction linear in rj. 
The modified self-similar solutions provide an excellent 



approximation to numerical solutions obtained for co ~ 3 
at large R, with 8 — > (and r\ 5= 0) for 3 < u> < u g . 

The second reason for revisiting the strong explosion 
problem is that it is of general methodological interest. It 
demonstrates transitions between (mathematically and 
physically) different types of solutions as the value of 
co changes, as illustrated in Figure [T] as the value of 
co increases above co = 3, a transition occurs between 
first-type s olutions (at co < 3) and se cond-type solutions 
(at co > 3; I Waxman fc Shvartsl IT993') ; as the value of co 
increases above co = ui c (l) ( w c ~ 8 for 4/3 < 7 < 5/3), 
the co < co c power-law solutions (R oc t 1 ^ 1- s \ 8 < 1) are 
replaced with exponential solutions (i? oc e*/ T , 8 = 1) at 
co = co c and with solutions diverging in finite time (R oc 
(-t) 1 ^*- 1 ), 8 > 1) at co > co c (jWaxman fc Shvartsll20Tfj|) . 
We show here that the strong explosion problem also 
exhibits a transition between two sub-types of second- 
type solutions at co = co g (7) . We argue in Section [4] 
that, based on the results presented in this paper, the 
definition of the two types of self-similar solutions should 
be somewhat modified and that the family of asymptotic 
second-type solutions should be expanded. 

Finally, we note that the propagation of shock waves 
in steep density gradients is of interest in a wide variety 
of astrophysical con texts (e.g. lOstriker fc McKeel 119881 : 
iKoo fc McKedfl990l and references therein), such a s su- 
pernova explosions (e.g. iMatzner fc McKeel I1999L and 
references therein). It is worth noting that self-similar 
solutions for shock propagation in power-law density pro- 
files, p oc r~", are useful for describin g shock propagation 
in more general density profiles (e.g. M atzner fc McKeel 
Il999tl0ren fc S ari 2009), as well as for the general study 
of shock wave stability ( e .g. iGo odman 1990; IChevalien 
119901: iKushnir et al . 2005; Sari et all 



2. SELF-SIMILAR SOLUTIONS IN THE "GAP" REGION 

As explained in Section [U we expect the strong explo- 
sion flow to approach a self-similar behavior, of the form 
given by Equations (JT]) and (JSJ) , as R diverges. Since 
for strong shocks the density just behind the shock wave 
is a constant factor, (7 + l)/(7 — 1), times the density 
just ahead of the shock, we must have e = — w, and 
we may choose B = K. With this normalization, the 
Rankine-Hugoniot relations at the shock front determine 
the bou ndary conditions for the s elf-similar solutions to 
be fe.g. lZel'dovich fc Raizerlll967f ) 



U(l) 



7 + 1 



c{1) = ^ G(1) 



7 + 1 



7-1 

The only parameter of the self-similar solution that re- 
mains to be determined is 8. 

2.1. ST and WS solutions 

The self-similar solution, given by Equations ([1]) 
and ([2]), depends on two independent dimensional con- 
stants, A and B = K. In the ST analysis, it is as- 
sumed that the second dimensional constant, in ad- 
dition to K, that determines the self-similar solution 
is E. In this case, dimensional considerations imply 
R oc (E^/K) 1 ^ 5 -^, i.e., 



(4) 
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For uj < 3 we have S < 0, i.e., decelerating blast waves. 
The flow properties are qualitatively different in the 
regimes uj < w vac = (7 — j)/('J+ 1) and uj vac < uj < 3 
(see Figure [2]). For u < u; vac , U tends to I/7 and C 
tends to infinity as £ tends to zero. For w vac < uj < 3, 
the self-similar solution contains an "evacuated" region: 
there exists some finite £; n > 0, such that the spatial re- 
gion < £ < £in is evacuated (p = 0). The self-similar 
solution describes the flow for £; n < £ < 1 and is matched 
to the evacuated region, £ < £; n , by a weak discontinuity, 
which lies at £ = £ ln . In this case, U tends to 1 and C 
tends to as £ tends to £j n . A detailed discus sion of the 
ST so lutions is given by iKorobeinikovi ()1991l ) and iBookl 
(fl99l . 

As explained in detail in ([Waxman fc Shvartsl fl993). 
the ST solutions are the correct asymptotic solutions 
only for u) < 3, for which 5 < 0. For larger values of 
uj the mass and energy contained in the self-similar so- 
lution are infinite, reflecting the fact that the initial gas 
mass at r > d diverges for d — > 0. I t was therefore 
suggested by IWaxman fc Shvartsl (|1993h that for uj > 3 
the asymptotic solution is given by a self-similar solution 
only over part of the (£, i?)-plane, bounded by £ = 1 and 
£ c (i?) < 1, and by a different solution at < £ < £, C {R)- 

Since such a solution includes a contact or a weak dis- 
continuity at £ C (R), £ c (i?) must coincide with a charac- 
teristic of the self-similar solution. For the self-similar 
flow, the characteristic lines 



C 

are given by 
C 
C±: 



dr _ 
dt ~ 

dlnR 
dln£± 

dhiR 



C± 



dr± 
~~dT 



u ± c 



(5) 



Ufa) - 1, 

= U&±) ± Cfa) - 1. 



(6) 



The directions in which the different characteristics prop- 
agate are illustrated in Figure [31 The physical interpreta- 
tion of the behavior illustrated in this figure is as follows. 
The flow just behind the shock is always subsonic: the 
shock-front point (U,C) — (f7(l),C(l)) lies above the 
"sonic line" U + C = l, i.e., [7(1) + C{1) > 1, which im- 
plies that C+ characteristics emerging from points just 
behind the shock always overtake it. C+ characteris- 
tics that do not overtake the shock exist only if the self- 
similar solution crosses the U + C = 1 line in the (U, C)- 
plane into the region where U+C < 1. Co characteristics, 
however, never overtake the shock and propagate away 
from it (in £ space). 

Requiring £ c to coincide with a C+ characteristic, that 
does not overtake the shock, implies therefore that the 
solution must cross the U + C = 1 line. Since A = 
along the sonic line, Equation (|A3|) imply that a physical 
solution must cross the sonic line at a singular point Ai = 
A2 = 0, as otherwise U(£) and/or C(£) are not single 
valued. This requirement determines the correct value 
of 5 for the uj > 3 asymptotic solutions. 

The self-similar solutions obtain ed in this way for 
uj > 3 wer e analyzed in detail in (jWaxman fc Shvartsl 
ll993ll20Toh . They describe accelerating blast waves, with 
S > 0. As £ — > the C(U) curves of these solutions ap- 
proach the singular point (U,C) = (1 — 6,0) for S < 1, 
and the singular point (U,C) = (0,0) for 6 > 1 (see 
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u < w v ac, S = Sst 

"vac < w < 3, S = S ST 

(jj g < uj < u) c , 5 = 6yps < 1 

lj c < u) : 5 = Sws > 1 

Sonic Line 

□ Strong shock point 

o Singular points 




Fig. 2. — Different types of the C(U) curves of the solutions of 
the strong explosion problem, obtained for d values in the regimes 

OJ < Uvac = (7 — 7)/ (7 + 1), Wvac < UJ < 3, LUg < UJ < Ulc, and UJc < 

uj. Also shown is the sonic line, U + C = 1. The square denotes 
the strong shock point, Equations (3), and the circles denote the 
singular points (U,C) = (0,0), (U,C) = (1 - 6,0), and (U,C) = 
(1,0). 
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Fig. 3. — Zoomed version of Figure[2] (using the same line types). 
Solid arrows indicate the propagation direction of C+ character- 
istics and dashed arrows indicate the propagation direction of Co 
characteristics. 



Figure [5]) . Analyzing the behavior of the solutions near 
these singular points, it was shown that although the 
mass and energy contained in the self-similar solution 
are infinite, the mass and energy contained within the 
region £ C (R) < £ < 1 approach finite values as R — > oo, 
for any C+ characteristic which satisfies £ C (R) —> as 
R — >• oo. Moreover, it was shown that £ c (i?)i? tx t as 
R — > oo, implying that the asymptotic flow within the 
region < r < £ C (R)R is described by the self-similar so- 
lution of expansion into vacuum. Finally, it was demon- 
strated by numerical calculations that the asymptotic be- 
havior described above is indeed approached for R/d ^> 1 
(IWaxman fc Shvartsl fl993l [20Toh . 

The self-simil ar solutions derived by 

IWaxman fc Shvartsl ()1993l ) exist only for lu > ^(7) > 3. 
Within the range 3 < uj < ui g (pf), there is no value of 
5 for which the C(U) curve crosses the sonic line at 
a singular point (oJ g ('y) > 3 is increasing with 7, with 
uj g = 3.26 for 7 = 5/3 and ui g — > 3 for 7 — > 1, see 
Figure [lj. Since the ST solutions provide the correct 
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asymptotic solutions only for uj < 3, the asymptotic 
behavior within the (narrow) range of 3 < uj < ui g (-^) is 
not described by cither of the two types of solutions. 

2.2. The asymptotic solutions for uj values within the 

"gap" 

As explained in the previous section, the divergence 
of the energy enclosed in the self-similar solutions for 
uj > 3 suggests that the asymptotic solution should be 
composed of a self-similar solution describing the flow 
at £ c (i?) < £ < 1, matched to a different solution at 
< £ < £c(-R)- In general, £ c ma y be a characteristic 
line o f any type. In the analysis of (jWaxman fc Shvartsl 
1993), it was assumed that £ c is a C + characteristic. 
This was mainly motivated by the fact that requiring 
the existence of a C+ characteristic that does not over- 
take the shock front as R — > 00 is equivalent to requiring 
that the solution passes through a sonic point, and it 
is commonly accepted that the similarity exponents of a 
second-type solution are determined by the requirement 
that the solution passes through such a singular point. 
We argue here that in order to determine the similarity 
exponents it is sufficient to require the existence of any 
characteristic that does not overtake the shock and for 
which the energy contained in the self-similar part of the 
flow, £ c (i?) < £ < 1, does not diverge as R — > 00. Ex- 
amining ui values below, within, and above the "gap" we 
show that for each value of u there is only one value of 
8 which yields a valid solution: a solution that does not 
cross the sonic line at a non-singular point and that either 
contains a finite energy or a characteristic line satisfying 
the conditions described above. 

The energy contained in the £ C (R) < £ < 1 region of a 
self-similar solution is 
rR 



E 3 (R): 



■ I 1 

drinr —pu 



£c{R)R 
-4ttR 3 



1 



7' 



1 



J+2S a2 



A'K{I k [UR)]+h[UR)]}, (7) 



with 



1 1 
4(0 = J dt'^G-U 3 , li(0 = J dt^G^—C 2 . 

The Ik and Zj terms give the kinetic and internal energy 
of the gas, respectively. In order for E S (R) not to diverge 
as R — > 00 we must have 

6<6 ST = (9) 

For S > Sst the energy contained in any £1 < £ < £2 
region of the self-similar solution diverges since j{ 3 - UJ + 2d 
diverges as R — > 00. For the WS solutions, 8 < Sst and 
7fc(£) diverges as £ — » in such a manner that the prod- 
uct i? 3_ " +2<5 /fc [£ C (-R)] tends to a constant as R diverges. 

Let us consider the C(U) curves, the behavior of char- 
acteristic lines, and the behavior of E S (R) for self-similar 
solutions with e = — uj, that satisfy the strong shock 
boundary conditions, Equations ([3]). We examine uj val- 
ues below the gap, w va c < w < 3, within the gap, and 
above the gap, uj g < to < uj c . Table [T] summarizes the 
relevant properties of the solutions obtained at the dif- 
ferent uj regions for different values of 8. The derivation 
of these properties is described below. 



First, we consider the properties of the C(U) curves by 
numerically investigating the solutions of Equation (|A2|) . 
We find that for 8 values smaller than a critical value, 
8*(w), the C(U) curves starting at the strong shock point 
cross the sonic line at a non-singular point. We find that 
8* = Sst for w vac < uj < 3, = for 3 < uj < ui g , 
and <5* = Sws f° r ui g < ui < ui c . As explained above, 
the self-similar solutions obtained for 8 < 8*(u>) are not 
physical. For ui g < ui and 8 = 8*(uS) = <5wSj the C{U) 
curves cross the sonic line at a singular point. Finally, for 
other values of 8, 8 > S*(uj) for 3 < uj < uj g and 8 > 8*(oj) 
for uj g < uj, the C(U) curves do not cross the sonic line 
and terminate at the singular point (U,C) = (1,0). 

Next, we examine the energy contained in the self- 
similar solution. For u) < 3, the only physical solution is 
that obtained for 8 — Sst- 8 must satisfy Equation ©, 
8 < 8st, in order for the energy in any £1 < C < £2 
part of the solution not to diverge, and the C(U) curve 
crosses the sonic line at a non-singular point for 8 < 8st- 
In order to examine the energy content of solutions in 
the (3 < uj < uj g ,8 > 0) and (uj g < uj,8 > <5ws) regions 
we need to analyze the behavior of the solutions near 
the singular point (U,C) — (1,0). The C(U) curves of 
the solutions in these regions of (uj, 8) lie above the sonic 
line, which implies that all C+ characteristics emerging 
behind the shock overtake the shock. Co and C_ charac- 
teristics, on the other hand, move along the C(U) curve 
toward the (U,C) = (1,0) singular point (see Figures [2] 
and El). 

Defining / = 1 — U, Equation (|A2j) is given, to leading 
orders in / and C, by 

j 1ti f f s f +[3-(n-2i) /7IC 2 for 5^0 



d\nC 



f 2 -(3-u/i)C 2 



(7-l)/M(7-lW(27)]C 2 



for 5 = 0. 

(10) 

Let us consider first the 8 > case. Regardless of 
whether / approaches faster or slower than C 2 , Equa- 
tion (|10|) implies that in this limit 



r dlnf 

llm ,, 7; 

/->o din C 



(11) 



where v is some constant. Assuming that / approaches 
zero slower than C 2 , i.e., that v < 2, leads to contra- 
dictions since Equation (fTOj) gives v = —2/(7— 1) < 0. 
Assuming that / approaches zero faster than C , i.e., 
that v > 2, Equation dTUJ> gives v = (67 - 2uj + 48)/ [(7 - 
1)uj + 2(5], for which v > 2 implies uj < 3. Thus, for 
uj > 3 and 8 > we must have v = 2, i.e., / oc C 2 , and 
Equation (fTUj) gives 



■> s 



-G 1 



(12) 



Using this result, the quadrature, Equation (|A3[) . gives 



/= 3(7-l) + 2* / X 

7 Uin 



(13) 



i.e., the singular point (U,C) = (1,0) is approached for 
finite £i n > 0. Using these results and Equation (|A4[) we 
find 

G OC J-(7"+2<5-3)/[3(7-l)+2<5]^ 

We may now determine the R dependence of the energy 
E S (R), contained within the £,o(R) < £ < 1 region of 
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Table 1. Determining S Based on the Properties of the Self-similar 
Solutions 



Properties of the self-similar solutions 


Wvac < O) < 3 


3 < uj < 0)9(7) 


"9(7) < w 


C(U) crosses the sonic line at a non-singular point 


S < <5sx 


<5 < 


<5 < <5ws 


C(U) crosses the sonic line at a singular point and 
E a (R) does not diverge as R — ¥ oo 






<5 = <5 W s 


C{U) terminates at (U,C) = (1,0) and 
E S (R) does not diverge as R — > oo 


S = <5 ST 


<5 = 




C(U) terminates at (U, C) = (1, 0) and 
E S (R) diverges as R — > oo 


5 > <5st 


s > 


5 > <5ws 



the self-similar solution, where £0 is a Co characteristic. 
Equation ([6]) may be solved, using Equation (fl"3|) . to give 



Using these results and Equation (|A4I) we find 



Sin 



_ R -[3( 7 -l)+25]/ 7 



(15) 



for the evolution of Co characteristics. Using Equa- 
tions (fTB")) and (O, we find that the kinetic energy in- 
tegral, lk(0 given in Equation ((SJ, diverges in the limit 
£ -> 6n as 

4(0 OC /(^)l-(-y«+2J-3)/[3(7-l)+2«] 

-7(o»-3)/[3(7-l)+2«] 

• (16) 



Using Equation (|T5|) we then find 

/felCo^oc^- 3 ), (17) 

which, using Equation ([7[l, implies that the energy di- 
verges (for 5 > 0) as E s oc i? 2<5 . 

This divergence of energy implies that the S — Sws > 
solutions are the only physical solutions for uj > u g . 
These solutions cross the sonic line at a singular point 
and terminate at (U, C) = (1 — 6, C = 0) with finite en- 
ergy E s , while the 5 > #ws > solutions terminate at 
(U, C) — (1, 0) with diverging energy E s . The divergence 
also implies that the 6 = solutions are the only solu- 
tions that may be physical solutions within the gap. Let 
us consider therefore the behavior of the 5 = solutions 
next. 

The solution of Equation (flO| for S = must also be 
of the form given by Equation (|11[) . Assuming / tends 
to slower than C, i.e., v < 1, leads to a contradiction 
since Equation (fT0|) gives ^ = 1/(7— 1) > 1. Therefore, 
v must satisfy v > 1. For v > 1, Equation (fT0|) gives 

67- 2w 



which satisfies v > 
Equation (fTUt gives 

f- 



(7-l)u/ 
1 for w < 67/(7 + !)• 



67 _ (7 -f 



c- 



27(2-7) 

The solution of the quadrature, Equation 



'In [ /- 

. Sin 



If" 

3(7- 



(18) 
For 1/ = 1, 

(19) 



(20) 



which implies 



In I ^ 



7+1 



3j), gives 
v > 1 ; 
i/=l, 



(21) 



Gck/- 



(7-1)^/(37-0;), v>\ 
((7 + l)w-6)/3( 7 -l), i/ = l, 



(22) 

which implies that ifc(C) diverges in the limit £ — > £i n as 
4(0 oc /(O 1 ^, which gives 



I k (to)<xR e ^-V = j?-3. (23) 

Numerical integration of Equations (IA2[) and (|A3I) shows 
that solutions starting at the strong shock point, Equa- 
tions approach the (U,C) = (1,0) singular point 
along a v > 1 curve. The asymptotic behavior of Ik is 
the same, as we find here, for both v = 1 and v > 1. 

Using Equation ([7]). we find that the kinetic energy 
part of -Eg approaches a finite, non-zero, constant as R 
diverges (it is straightforward to verify that the internal 
energy part of E a vanishes in the limit R — > 00). The 
fact that the kinetic energy approaches a constant also 
implies that the mass contained within £o(R) < £ < 1 
approaches a finite non-zero constant, since the velocity 
of each fluid element, £o-R, approaches a constant, ^ ln R. 
Thus, the 5 = solutions satisfy the requirement for the 
existence of a characteristic £ C (R) that does not overtake 
the shock, and for which the energy contained in the self- 
similar part of the flow, £ c (i?) < £ < 1, does not diverge 
as R — ^ 00. 

As pointed o ut in Section [U it was suggested by 
iGruzinovl (|2003l) that the asymptotic solutions in the gap 
are the 5 = solutions. The justification given there 
is based on the argument that the non-self-similar part 
of the asymptotic flow, which must exist since the self- 
similar solution contains infinite energy, acts as an infi- 
nite mass piston, which must move at a constant speed 
and may support the flow ahead of it. As we show here, 
the mass and energy contained in the self-similar and 
non-self-similar parts of the flow are both finite and may 
be comparable (this is also the case for the WS solutions 
obtained for u > uj „). Thus, the va lidity of the heuristic 
argument given in "(Gruzinov 2003) is not obvious. 



3. SLOW CONVERGENCE TO SELF-SIMILARITY: 
NUMERICAL RESULTS AND MODIFIED SELF-SIMILAR 
SOLUTIONS 

3.1. Comparison to numerical solutions 

We present in this section a comparison between the 
5 = self-similar solutions and numerical solutions of 
the flow equations obtained for u> values within the gap. 
We have numerically calculated the propagation of a 
strong spherical shock wave through an ideal gas using a 
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w = 3.1, 7 = 5/3 



0.75 



Strong shock point 
Sonic Line 

Numerical R/d = 9 • 10 4 
Self-similar 5 = 
Self-similar 5 = -0.01865 




0.85 0.9 

U 



Fig. 4. — C(f7) curve for uj = 3.1, 7 = 5/3. Shown are the 
numerical solution at R/d = 9 X 10 4 and two self-similar solutions, 
corresponding to <5 = and 8 = —0.01865. Note that the C(U) 
curve of the <5 = —0.01865 self-similar solution crosses the sonic 
line (at a non-singular point). This is not visible in the figure since 
the crossing takes place very close to ((/, C) = (1, 0). 



-0.018r 




-0.025 
-0.026 
-0.027 



x 10 



Fig. 5. — <5(iJ) determined from a numerical simulation, using 
the zero slope method (see the text), for an explosion with w = 
3.1, 7 = 5/3. The error bars are an estimate of the accuracy of 
the determination of 8 using this method. The accuracy is better 
for larger values of R/d, where the flow behind the shock is better 
resolved. 



Lagrangian sche me with total energy conservation (e.g. 
ICaram ana 1998). Shock waves are described in our cal- 
culations using a von Neumann artificial viscosity, im- 
plemented as a pressure term in all cells with negative 
difference between the nodes' velocities, Au < 0, 

q = Au (x q Au — xic) , (24) 

with x q — 4 and x\ = 0.1. The initial conditions used are 
zero velocity everywhere, constant density and pressure 
at r < d, relatively small pressure at r > d (see below), 
and a density profile proportional to r~ u at r > d. The 
initial mesh spacing was uniform, Ar — d, i.e., the pres- 
sure was high only in the innermost cell. We chose the 
density of the innermost cell such that its mass is 10 times 
higher than its neighbor's, and we chose its pressure to 
be 10 6 times higher than its neighbor's. The pressure 
at the rest of the cells was chosen such that the outgo- 
ing shock wave is always strong. The mesh included 10 5 
cells, i.e., we were able to calculate explosions up to a 
radius of R/d ~ 10 5 . 

We would like to examine the behavior of S(R) = 
dhiR/dlxiR = RR/R 2 . Since 8 depends on R and the 
numerically determined value of R is noisy due to the 
finite resolution, a derivation of 6 by a direct differen- 
tiation of R is not accurate enough for our study. In 
order to overcome this problem, we derive 8 from the 
numerical spatial profiles: we choose S as the value for 
which the difference between the numerical profiles and 
the self-similar profiles, determined by Equations (|A2I) - 
(|A4I) . has a zero slope at the shock front (this ensures 
that the acceleration of the shock wave, which is deter- 
mined by the spatial profiles in the vicinity of the shock 
front, is the same for the self-similar and numerical so- 
lutions). In what follows, we give some details regarding 
this method for the determination of 8, which we term 
the "zero slope method" . 

We examine the difference /(£) between the self-similar 
profiles (U, C, and P), obtained for a chosen value 
of 8, and the profiles obtained in the numerical simu- 
lations at some R/d. For the comparison of the self- 



similar and numerical profiles we consider the radial 
range 0.999 > £ = r/R > 0.9995. The upper limit of 
this range is set by requiring that the oscillations behind 
the shock, caused by the artificial viscosity, are damped 
considerably, while the lower limit is chosen to ensure 
that deviations from the self-similar solutions are small 
(see Section |3.2|) . In order to determine whether or not 
the difference /(£) is consistent with a zero slope, i.e., 
with /' = 0, it is insufficient to use a simple linear fit for 
/(£), since oscillations in the numerical profile caused by 
the artificial viscosity are not random and cannot be ne- 
glected. We define therefore / = (/ — and 
£ = (£ — where /i and a stand for mean 

and standard deviation over all numerical cells in the 
range, and examine the number of points for which £ > 1 
and / > 0(< 0), denoted 7V ++ (_), and the number of 
points for which £ < 1 and / > 0(< 0), denoted by 

N U—)- In the absence of numerical inaccuracies, the 

ratios r + r_) — N + ^_- j+ /N + ^_- ) _ should all equal unity for 
/' = 0. In order to allow for numerical inaccuracies, we 
consider /' to be consistent with for 0.1 < r+(_) < 10, 
and determine the range of allowed values of 8 as the 
range for which 0.1 < ?"+(-) < 10. Since »"+(_) (or 
l/r + (_)) grow rapidly as 8 is modified, the range of al- 
lowed values of 8 is not sensitive to the exact choice of 
the allowed range of r+r_). 

In Figure |4j we compare the numerical C(U) curve 
obtained at R/d — 9 x 10 4 , for an explosion with 
u> = 3.1, 7 = 5/3, with the self-similar curve obtained 
for the value of 8 determined by the zero slope method, 
8 = —0.01865. The numerical curve and the self-similar 
one are very close near the shock front and show small 
discrepancy far from the shock. Similar results are ob- 
tained for other values of R/d. 

S(R) determined by the method describe above is 
shown in Figure[5]for an explosion with 7 = 5/3, u) = 3.1. 
The error bars are an estimate of the accuracy of the 
determination of 8 using this method. It is apparent 
that the convergence of the inferred value of 8 is very 
slow, and that it has not converged for a very large value 
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of R/d, R/d » 10 5 . It is difficult to determine, based 
on the simulation, whether or not 8 approaches for 
R/d — > oo: the rate of change of 5 obtained at R/d ~ 10 5 , 
dd/dln(R) «3x 10~ 3 implies that an increase in R/d 
by a factor of 10 will modify the inferred value of S from 
S w -0.019 to only give S « -0.014. 

A comment is in order regarding the convergence of our 
numerical solutions. We have checked convergence by 
examining the modification of the inferred value of 5 at 
fixed R/d, obtained when the number of cells within the 
region of high initial pressure is increased, i.e., choosing 
Ar/d < 1 (and keeping a uniform initial grid spacing). 
The results presented in this section for S at large R/d 
using Ar/d = 1 are converged to a few percent. For 
example, for the 7 = 5/3, u = 3.1 case, increasing the 
number of cells by a factor of 2, i.e., using Ar/d = 1/2, 
changes the inferred values of S by less than a 5%. 

Slow convergence of the numerical solutions to an 
asymptotic self-similar behavior is obtained also for other 
values of 7, lo within or near the gap (see also Section|4]). 
This result motivates us to explore in Section 13.21 mod- 
ified self-similar solutions that describe the approach of 
the flow to self-similarity. 



■ 3.1, 7 = 5/3 



3.2. Modified self-similar solutions 
To quantitatively consider the approach to 



self- 



similarity, we examine solutions of the hydrodynamic 
equations of the form 

u(r,t)=Rt[U(Z,6(t)) + r,U 1 (Z,5(t))], 

c(r,t) = M[C^,d(t)) + r ) C 1 (^8(t))], 

p(r,t) = BR s [G(£,5(t)) + r 1 G 1 (Z,5(t))], 

p(r, t) = BR £ R 2 [P(Z, S(t)) + ijPifc S(t))}. (25) 

Here, 5(t) = dlnR/dhxR, r, = d5/dhxR, and F(£,5{t)) 
stand for the self-similar solution F(£) obtained for the 
instantaneous value of 6, S(t). The flow fields are of the 
general form 

f(r,t) = R a RP[F(£,6(t)) +vFi&8{t))] 
For this form, 

=R a - 1 &(F' + V Fl), 

where ' = d/d£, and, using 5 — f]R/R, 
Of 



(26) 



(27) 



-± = a R^R^m, 6(t)) + nF, (£, S(t))} 



■f3R a RP~ 1 R[F(£, S(t)) + r]Fi(£, 6(t))} 
.R"Re[(F' + r,F{(Z,5(t)))(-^R) 
dF 



85 



dF x ■ 
V^S + ijFx] 



R^RP+^F + r)F 1 ){a 
,8F dF x dlnr) 

Restricting to solutions with dlnrj/dhiR 



p8) 
-Fi)]. 



Z(F' 



rjF{) 

(28) 
and 



neglecting 77 terms, assuming 77 = dS/d\nR <S 1, 
we obtain a set of ordinary differential equations for 
Ui , Cx, Gi, which may be written as 

Ux\' /Ux\ fl / U 

(29) 




1 Sonic Line 

Numerical, R/d = 9 ■ It) 4 
-6 = -0.01865, 17 = 
-8 = -0.01865, 7] = 2.13-10- 




0.85 0.9 

u = 2.8, 7 = 5/3 



' ' ' Sonic Line 
o Numerical, R/d = 9 ■ 10 4 
- -S = -0.1029, n = 

S = -0.1029, rj = 0.6- 1Q- 




FlG. 6. — C(U) curves for explosions with u = 3.1,7 = 5/3 and 
w = 2.8, 7 = 5/3. Shown are the numerical solutions, self-similar 
solutions (7) = 0) and modified self-similar solutions (7] 7^ 0). The 
difference between the r\ = and rj 7^ solutions is d ifficult to 
identify in this plot. It is clearly shown in Figures [611101 
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Fig. 7. — Numerical and modified self-similar spatial velocity pro- 
files of the calculations presented in Figure [6] Top panels present 
a comparison of the numerical profiles with the self-similar pro- 
files o btai n ed f or the appropriate value of 5 (solutions of Equa- 
tions l)A2j l- llA4| l'). Bottom panels present a comparison of the dif- 
ference between the two, i.e., of the numerical deviation from the 
self-similar solution, with the deviation pre dict ed by the modified 
self-similar solutions, defined by Equations i ]2-">[ ) and determined by 
Equation J29j. 



The matrices A, B, and C are given in Appendix [Bj The 
equations for U, C, G are the same self-similar equations 
as before, given in Appendix [X] The boundary condi- 
tions at the shock are set by the Rankine-Hugoniot re- 
lations, which imply 

UMt), 1) = P x {&(t\ 1) - Gi (<*(*), 1) = 0. (30) 

Note that n does not appear in the equations describ- 
ing the modified solutions, Equation (|29|) . and cannot 
therefore be determined by these equations. 

We compare in Figures IMTOl the modified self-similar 
solutions to the results of numerical calculations, for 
oj = 3.1, 7 = 5/3 and u = 2.8, 7 = 5/3. We have chosen 
an example with u) = 2.8 < 3 below the gap, in order to 
emphasize that slow convergence to self-similarity is not 
unique to the gap region. The C(U) curves are compared 
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Fig. 8. — Same as Figure [7] but for the spatial pressure profiles. 
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Fig. 10. — Same as Figure[7] but for the spatial density profiles. 
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Fig. 9. — Same as Figure [7] but for the spatial sound velocity 
profiles. 



in Figure [BJ and the spatial profiles are compared in Fig- 
ures [THUS The spatial profiles U\, C\, Pi, G\ of the nu- 
merical solutions were obtained by subtracting from the 
numerical profiles the self-similar profiles U, C, P, G cor- 
responding to the appropriate value of 5(R). The value 
of rj was determined by comparing the numerical results 
for Ui, G\, Pi, Gi with the solutions of Equation (|29[) 
(note that the solutions of Equation (f29| are indepen- 
dent of rj, which just determines the normalization of 
the deviation from the self-similar solution, as given by 
Equations (|25]l ). The agreement between the numerical 
solutions and the modified self-similar solutions, demon- 
strated in the plots of Figures [BHTU1 suggests that the 
modified self-similar solutions provide an approximate 
description of the approach to self-similarity. 

4. SUMMARY AND DISCUSSION 

We have shown that the self-similar solutions describ- 
ing the asymptotic flow of the strong explosion problem 
for lj values within the gap, 3 < lj < Lj g (^), are the 
5 = 0, R = constant, self-similar solutions. For lj > 3, 
the energy in the self-similar solutions is infinite, imply- 
ing that the self-similar solution may describe the flow 
only in part of the (£, i?)-plane. This suggests that for 



lj > 3 the asymptotic solution should be composed of a 
self-similar solution describing the flow at £, C (R) < £, < 1, 
matched along some characteristic line £, C (R) to a differ- 
ent solution at < £ < £ C (R), such that the energy con- 
tained in the self-similar par t of the flow, £ C (R) < £ < 1, 
does not diverge as R — > oo (IWaxman fc Shvartill993D . 

In the analysis of (|Waxman fc Shvartsl I1993T T it was 
assumed that £ c is a C+ characteristic. This was mainly 
motivated by the fact that requiring the existence of a 
C+ characteristic, that does not overtake the shock front 
as R — > oo, is equivalent to requiring that the solution 
passes through a sonic point, and it is commonly ac- 
cepted that the similarity exponents of a second-type so- 
lution arc determined by the requirement that the solu- 
tion passes through such a singular point. We showed 
here that in order to determine the similarity exponents 
it is sufficient to require the existence of any character- 
istic, that does not overtake the shock, and for which 
the energy contained in the self-similar part of the flow, 
£ c (i?) < £ < 1, does not diverge as R — > oo. 

Examining lj values below, within, and above the 
"gap" , we showed that for each value of w there is only 
one value of S, S = <5*(w), which yields a valid physical 
self-similar solution (see Table [T]). For S < <5* the C(U) 
curve determined by Equation (|A2[) crosses the sonic line 
U + C = 1 at a non-singular point (yiel ding a non- 
single- valued solution, see Equations (|A3[) and (|A6p ). 
For S > <5*, the self-similar solution energy diverges, in 
the sense that the energy contained in £, C (R) < £, < 1 
diverges for any choice of a characteristic, £ c (i?), that 
does not overtake the shock wave. The physical solutions 
are the ST solutions, 6 * = Sst = {uj — 3)/2 for a; < 3 , 
the solutions derived in (jWaxman fc Shvartsl fT993. 2010), 
S* = <5\vs < <5st, for lj > uj g , and the 5* = solutions 
for 3 < u) < ojg. The C(U) curves of the 3 < to < u g 
solutions do not cross the sonic line and £ C {R) must be 
chosen as a Co characteristic for these solutions. 

In Section [3~T1 we compared the asymptotic, R/d » 1, 
behavior of numerical solutions of the hydrodynamic 
equations, Equations (|Al|) . to that expected based on 
the (5 = self-similar solutions. We find that while the 
flow approaches a self-similar behavior with \S\ <§; 1, the 
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Fig. 11. — Self-similar exponent 6 as a function of u> for 7 = 5/3 
(solid line), compared with the value of 8 inferred f rom numerical 
simulations (using the method described in Section 13.21 1 at R/d = 
9 X 10 4 (dashed). 



convergence to self-similarity is very slow for oj ~ 3 (e.g., 
Figure [5]) . It should be noted that convergence to self- 
similarity is slow for any value of u ~ 3, both within and 
below the gap, as demonstrated in Figure [TT] Hence, 
it is difficult to check using numerical solutions whether 
for uj values within the gap the flow indeed approaches 
a S = self-similar behavior as R — > oo. We showed in 
Section 13.21 that in this case the flow may be described 
by a modified self-similar solution, d In R/d In R = 8 with 
slowly varying S(R), rj = dS/d\nR <C 1. In these solu- 
tions, the spatial profiles are given by a sum of the self- 
similar solution corresponding to the instantaneous value 



of S and a self-similar correction linear in 77, see Equa- 
tions (|25|) . The equations describing the self-similar cor- 
rections are given in Equation (|2"9")l . We have shown that 
for u ~ 3 the modified self-similar solutions provide an 
approximate description of the flow at large R, see Fig- 
ures EHH3 with s -> ( and 7^ 0) for 3 < w < These 
results support the conclusion that the flow approaches 
the S — self-similar solutions as R diverges. 

Based on the analysis presented here, we suggest that 
the definition of first- and second-type similarity solu- 
tions should be somewhat modified, and that the family 
of second-type solutions should be expanded. Solutions 
of the first-type may be defined as solutions that are 
valid over the entire (r, <)-plane (or the part of which 
where the flow takes place). Such solutions must satisfy 
the global conservations laws (of mass, momentum, and 
energy) , and hence the values of the similarity exponents 
of such solutions may be determined by dimensional con- 
siderations. Solutions of the second-type may be defined 
as solutions, which are valid only in part of the region 
in the (r, i)-plane over which the flow takes place. Such 
solutions should be required to allow the existence of 
a characteristic line, £ C (R), along which the self-similar 
solution is matched to another solution, and to comply 
with the global conservation laws within the region of the 
(r, i)-plane described by the self-similar solution. 
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APPENDIX 



THE EQUATIONS DESCRIBING SELF-SIMILAR FLOWS 
The equations describing adiabatic one-dimensional flow of an ideal gas are (e.g. lLandau fc Lifshitzl ll959) 

(d t + ud r ) In P + r-( n - 1 >«9 r (r"-y) - 0, 
(d t + ud r )u + p- 1 d r { 1 - 1 P c 2 ) = Q 1 
(dt+ud r )(c 2 p 1 -<)=0, 



(AI) 



where n = 1, 2, 3 for planar, cylindrical, and spherical symmetry, respectively. 

Substituting Equations ([H and (El) in the hydrodynamic Equa tions (|A1|) . the partial differential equations are 
replaced with (Zcl'dovich & Raizer 1967; Waxman & Shvarts 2010) a single ordinary differential equation, 



dU _ Ai(£/,C) 
dC ~ A 2 {U,C)' 



and one quadrature 



G is given implicitly by 



dln£ MU,C) dln£ A(U,C) 

or 



dU Ai (U,C) " dC A 2 (U,C)' 

(iCy 2{n+c) \\ - J7|A G (7-D(n+ e )+AgnA = const 



(A2) 



(A3) 



(A4) 



with 



A = -(7-l)e + 2& 



(A5) 
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The functions A, Ai, and A2 are 

A = G 2 - {l-U) 2 , 
/± x = U{l-U){l-U -5)-C~ [ nl' 
A 2 = C{(l-U)(l-U -6) 

' ~~ l -U[{n-l){l-U)+5]-C 2 



2 

2(J-(7-l)c c 2 



}■ (A6) 



2 7 1-U J 

THE EQUATIONS DESCRIBING MODIFIED SELF-SIMILAR SOLUTIONS 
The matrices A, B, and C of Equation (f2T)f are given, for spherical symmetry, by 

/ -£G \ 

A = eG(l -U) -1 (Bl) 

V - 7 £P(1 - U) £G(1 -U)J 

( 3G + £G' s + W + ZU' \ 

B= I £G(5 + 2U -1)+£ 2 U'G £U(5 + U -1)+(, 2 U'{U -1) (B2) 

V £GP'-^PG' P{e{l- 1 )+28)+£P'{U - 1) G(e(l - 7) + 25) - £G' 7 (?7 - 1) / 

and 

/ 1 \ 

C = £G . (B3) 
V -7P G / 
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